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Abstract 



o 

I Using the Hartree-Fock approximation, we calculate the energy of different 



Wigner crystal states for the two-dimensional electron gas of a double quan- 



tum well system in a strong magnetic field. Our calculation takes interlayer 

CO - 

hopping as well as an in-plane magnetic field into consideration. The ground 



state at small layer separations is a one-component triangular lattice Wigner 
state. As the layer separation is increased, the ground state first undergoes a 
transition to two stacked square lattices, and then undergoes another transi- 



tion at an even larger layer separation to a two-component triangular lattice. 



The range of the layer separation at which the two-component square lattice 
occurs as the ground state shrinks, and eventually disappears, as the inter- 

X" 

$— i ' layer hopping is increased. An in-plane magnetic field induces another phase 

transition from a commensurate to a incommensurate state, similar to that of 
v = 1 quantum Hall state observed recently. We calculate the critical value of 
the in-plane field of the transition and find that the anisotropy of the Wigner 
state, i.e.,, the relative orientation of the crystal and the in-plane magnetic 
field, has a negligible effect on the critical value for low filling fractions. The 
effect of this anisotropy on the low-lying phonon energy is discussed. A novel 
exerimental geometry is proposed in which the parallel magnetic field is used 
to enhance the orientational correlations in the ground state when the crystal 
is subject to a random potential. 
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I. INTRODUCTION 



An electron gas is expected to condense into a Wigner crystal [|T|(WC) below some 
critical density. This condensation occurs when the Coulomb energy, which tends to localize 
electrons into individual lattice sites to keep them as far apart as possible from each other, 
dominates over the kinetic energy, which favors a smooth variation of electron density In 
the absence of a magnetic field, the kinetic energy of a two-dimensional electron system 
(2DES) scales like K = h 2 /m*a 2 , while the Coulomb energy scales like V = e 2 /ea, where 
a is the mean inter-electron distance and e is the dielectric constant of the host material, 
and m* and e are the electron mass and charge respectively. The relevant parameter is the 
ratio r s = V/K = a/as, where clb = h 2 e/m*e 2 is the Bohr radius. Monte Carlo simulation 
predicts that a 2DES crystallizes for r s > 35. When a strong magnetic field is applied 
perpendicular to the 2DES, the situation is changed qualitatively, as the kinetic energy is 
quenched into discreet Landau levels, and the zero-point fluctuations in the lowest Landau 
level are confined within a magnetic length l Q = (hc/eB±) 1 ^ 2 , where c is the speed of light 
and B± is the applied magnetic field. Once l a is sufficiently small compared to the typical 
inter-particle distance a, crystallization occurs. The ratio l a /a can be characterized by the 
Landau level filling factor v = 2irl 2 n, where n is the density of the 2DES. Crystallization 
will occur for sufficiently small v for any given density. Theoretical estimates || put the 
critical filling factor of crystallization at about v c ~ 1/6. In recent years, especially after the 
observation of the reentrant insulating phase around the v — 1/5 quantum Hall state [|J, 
there has been an increasing interest in the study of the WC states in a 2DES in a strong 
magnetic field. Many experimental results PHlOfl are found to be in some ways consistent 
with the assumption of a pinned WC as the ground state. 

Recent advances in material-growth technology allows the fabrication of high quality 
double quantum well system (DQWS), in which two interacting 2DES are separated by a 
distance comparable to the mean inter-particle distance within the 2DES. This introduces 
a new degree of freedom associated with the third direction. Electron-electron interactions 
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between the layers have been known to lead to new quantum Hall states ||11|| . It also leads 
to increased stability of the WC state or other charge density wave states ||12j| . Recent 
experiments on two layer systems in wide quantum wells below filling factor v = 1/2 have 



shown insulating behavior |13| similar to that seen in single layer systems below v = 1/5. 
Because Coulomb interactions can lead to mixing of the electronic states of the two wells, 
more complicated structures of WC states become possible in DQWS. The goal of this paper 
is to investigate the evolution of the ground state among different possible WC configura- 
tions, as the parameters of the DQWS are changed, at a small Landau level filling factor 
where a WC is expected to be the ground state. This is accomplished by calculating, in 
the Hartree-Fock approximation, the ground state energy of the different WC states. Our 



method of calculation is based on the numerical technique developed in Ref. fT4 |, which is 



valid in the strong-field limit. We take into consideration inter-layer hopping as well as an 
in-plane magnetic field. We ignore the finite thickness of the quantum wells and treat the 
electron gas as ideally two-dimensional. This zero-thickness approximation is an important 
simplification in the case of a tilted field, since the effect of the in-plane magnetic field can 
then be included by adding a phase factor to the wavefunction of the electrons in one of the 
wells. 

At small layer separations, where the electron-electron correlations between the wells are 
almost as important as the correlations within each individual well, electrons will occupy 
the symmetric state of a DQWS to minimize Coulomb energy and to take advantage of 
hopping energy. In this situation, the DQWS behaves essentially like a single layer system. 
The WC at small filling factor is therefore a triangular lattice. As the layer separation 
increases, the 2DES will seek a state where intralayer correlations are favored over the 
interlayer correlation. This leads to a transition to a truly two-component system when the 
intralayer Coulomb energy is more important than the hopping and interlayer interaction 
energy combined. Under this condition, each individual well forms a WC separately with a 
lattice constant which is appropriate for the electron density in its own layer. The two lattices 
are shifted relative to one another to lower the interlayer static Coulomb energy (Hartree 



energy). We find that, in the absence of interlayer hopping, the one-component WC is first 
transformed to a two-component square lattice WC state, and then transformed to a two- 
component triangular lattice, as the layer separation is increased. The existence of a square 
lattice WC at intermediate separations can be understood by noticing that the square lattice 
configurations have lower interlayer interaction energy than that of the triangular lattices. 
This gain may exceed the difference in intralayer Coulomb energies between a square lattice 



and a triangular lattice, which is known to be small ||15|| . The range of the layer separation 
at which the two-component square lattice exists shrinks, and eventually disappears, as the 
interlayer hopping is increased, as a result of the expansion of the one-component triangular 



phase. Our result is consistent with that of Chan and MacDonald ||16|| , where they treat 
electrons classically. 

A DQWS in a tilted magnetic field has been studied recently. A new phase transition, 
driven by the in-plane component of the field, was observed for v = 1 quantum Hall state 
T7J, and was explained using an easy-plane itinerant quantum ferromagnetism description 
Tj|. This phase transition happens as a result of the competition between the hopping 
energy and Coulomb interaction energy. The in-plane magnetic field twists the interlayer 
phase coherence of the wavefunctions of electrons in symmetric state. The result is an 
increase in interlayer Coulomb energy. At small in-plane field, the increase in interlayer 
Coulomb energy is small compared to the hopping energy, so the electrons will stay in the 
state dictated by the tunneling part of the Hamiltonian in order to take advantage of the 
hopping energy. As the in-plane field becomes stronger, the hopping energy is reduced, 
while the cost in Coulomb energy continues to rise. When the in-plane field is raised beyond 
a critical value, the DQWS system will undergo a transition to a state in which electrons 
give up the hopping energy in order to restore the interlayer correlations. The same physics 
occurs for either a v — 1 quantum Hall state or a one-component WC states, since for 
both cases there is strong interlayer coherence before the application of an in-plane field. 
Since a WC breaks the rotational symmetry of the system, the phase transition should, in 
principle, depend quantitatively on the angle between the direction of the in-plane field and 



the crystal axes. We have calculated the critical value of the in-plane field and found that 
this dependence is very weak, practically unobservable at small filling factors. 

Nevertheless, a parallel magnetic field can have important consequences, particularly 
for the long-wavelength physics of the system. The broken orientational symmetry of the 
groundstate due to a parallel magnetic field implies that there is a restoring force on the 
crystal if it is misoriented with the field. This effect can be very important if the system 
is subject to a weak random potential, which destroys long-range orientational order in 



this system |T{|. We propose a novel experimental geometry, in which the 2DEG is cooled 
through its freezing transition in the presence of a parallel magnetic field. The force that 
tends to align the crystal axes perpendicular to the parallel field should enhance the orienta- 



tional order in the groundstate. It has been shown recently [19] that the depinning electric 



field [H of a disordered WC is sensitive to the orientational correlations in the system, and 
tends to be reduced as orientational order increases. Thus an increase in orientational order 
can in principle be detected experimentally through the transport properties of the system. 

This article is organized as follows. In Sec. 2, we describe the Hartree-Fock approxima- 
tion for the interacting 2DES in a DQWS in a strong magnetic field. In Sec. 3, we present 
and discuss the numerical results. Sec. 4 discusses the effect of a parallel magnetic field on 
some of the crystal properties. A brief summary in Sec. 5 concludes this article. 

II. HARTREE-FOCK APPROXIMATION 

In the absence of an in-plane magnetic field, the Hartree-Fock approximation for a DQWS 



has been clearly presented in detail in Ref. [14]. We need only to extend it to the case of 



a tilted magnetic field. In a Hartree-Fock approximation, one treats the interacting Hamil- 
tonian as that of free electrons in the mean-field potential determined by a given electronic 
state and then self-consistently solves for the state. In our present case, the electronic 
state being sought is a WC. The characteristic of a WC is a periodically modulated charge 
density. We therefore choose the Fourier transformed electron density at corresponding re- 
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ciprocal lattice vectors (RLV) as the order parameters of the WC states. Following this 
idea, we define Py(q) = / d 2 re~ tcl ' r ipj(r)ipj(r), where i,j = 1,2, labeling the two layers, ^ ( 
^■(r)) is electron creation (annihilation) operator. One obtains, in the lowest Landau level 
with the Landau gauge, 

p*M) = \ E e-^ a+ ^%, a+qy clc jP , (1) 

9 a/3 

where a and (5 are the single-particle states of the lowest Landau level. In the above 
expression, g = Vl/2itI 2 is the Landau level degeneracy, where Vl is the area of the DQWS. 
The order parameters have, by definition, the property that < Py(q) >=< Pji(— q) >*• We 
restrict ourselve to seek only the WC states where the charge distribution in one layer is 
the same as that of the other layer, except that one is rigidly shifted by a displacement of 
a relative to the other, i.e., < P22(q) >= e iqa < pn(q) >. This leaves us with only two 
independent sets of the order parameters to be obtained. The order parameters are related 
to Green's function as 

<p^(q) >=G J - i (q,r = 0-), (2) 

where G^(q, r) = / d 2 re _iq ' r Gy(r, r, r) and Gy(r, r',r) = - < T T ipi(r,T)ipj(r') >. In the 
lowest Landau level, the Green's functions become 

G y (q, r) = - ]T e-^ a+ ^ 2 5 a ^ +qy G tJ (a, (3, r), (3) 

9 a f3 

with 

Gij(a, P,t) = -< T T C ia {T)C)p > . (4) 

The Green's functions are to be obtained by self-consistently solving the equation of 
motion 

^G ij (a,P,r) + 5(r)5 ij 5 a , /3 + < T T [H,C ia {r)]C} p >= 0. (5) 

The next step is to approximate the equations of motion in the Hartree-Fock approach. 
We will treat a tilted magnetic field from the beginning. The result for a perpendicular 
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field is recovered by simply setting the in-plane component of the field B» = 0. Under the 
zero well-thickness approximation, the only effect of the in-plane magnetic field is to add a 
phase factor to the electronic wave function of one of the layers. Let k B = d/l 2 , with d the 
separation between the wells, and l\\ = (hc/eB^) 1 ^ 2 . If one chooses the direction of B\\ as 
the x-axis, the single-particle eigenstates become 

lQ ( r )= 1 e iay e -(x-al*r/2l2 



2a (r)= e ifc ^0 la (r). (6) 

The Hamiltonian of the DQWS contains kinetic energy, (which is a trivial constant for 
the lowest Landau level,) interlayer hopping energy, and the Coulomb interaction energy. 

H= -J2 nCj a C ia - te- k2 ^Y,(CtC2a- kB + Cl_ kB C la ) 

ia a 

+ ^EE V tJ (a,a\[3,P')Clc}fy,C ta ,, (7) 

where \x is chemical potential to be fixed at the end of calculation for a given electron density, 
and Vij(a, a', /3, (3') is the matrix element of Coulomb potential. The hopping parameter t 
is suppressed in the presence of an in-plane field by a factor e~ k B l °^ y which comes from the 
matrix element < la|t|2a' >= te~ k B l °/ A 5 atQ ' + k B - Physically, this means that the electrons 
tunnel along the direction of the total magnetic field. Performing the Hartree-Fock pairing 
C\CIC 2 ,C V ^< C\C V > C\C 2 ,- < C\C 2 , > C\C V on Eq. (0), one obtains 

e 2 

H = -g2_,iipu{ti) -fltf[P2i(k B ) +pi 2 (-k B )] +£-r2ZzZ* 7 y( c l)^( c l)> ( 8 ) 

i a o q ij 

where k B = k B y and 

Mq)= [K(q) - H(q)] < Pn(-q) > +K(q) < p 22 (-q) >, 

U 12 (q)= -V d (q) <p 21 (-q) >, (9) 

with U21 and U22 obtained by interchanging the indices 1 and 2. In the above expressions, V a , 
Vb, V c , and Vd are the direct and exchange terms for the interlayer and intralayer Coulomb 
interactions 



V a (q) = —e-^\l - 5 Qi0 ) 
V c (q) = e-« d V a (q) 

V d (q) = r dxJ (xql o )e' x2 ^- xd / 1 ", (10) 
Jo 

where Jo and Jo are the zero-th order Bessel function and modified Bessel function, respec- 
tively. It is worthwhile to notice that since Coulomb scattering is not expected to move 
electrons from one well to the other, the matrix element Vij(a, a', f3, f3') is unchanged by 
the presence of an in-plane field, i.e., Eq. ([|) and Eq. fllPf) are exactly the same as the ex- 
pressions for B\\ = 0. This means that the commutator [H c ,Ci Q ], where H c is the Coulomb 
interactions part of the Hartree-Fock Hamiltonian of Eq. is unchanged by the presence 
of the in-plane field. Letting Aij(a, (3) = — < T T [H c ,Ci a ](r)Cjp >, we have 

Aa(<h r) = ~ £ U ik (p - q)e^ A P/ 2 G fcj (p, r), (11) 
eL ° k P 

where q A p = q x p y — q y p x . The effect of the in-plane field is contained explicitly in H t , the 
hopping term in Eq. (§). Denoting Fij(a,j3) = — < T T [H t , Cial^Cjp >, one obtains 

F n (q, r) = te- k B l2j4 e- ll ^ kB/2 G 21 (q - k B , r) 
F 12 (q + k B , r) = te-^ / °/ 4 e-^^ fcfl / 2 G 22 (q, r) 
F 21 (q - k B , r) = te-W^'WGniq, r) 

F 22 (q, r) = te- k2 B l2/i e il ^ kB/2 G 12 (q + k B , r) (12) 

One can see, from the above expressions, that in the presence of an in-plane field the 
hopping Hamiltonian attempts to shift the positions of the non-vanishing interlayer order 
parameters < pi 2 > and < p 2 i > in reciprocal space from q = G to q = G±k B . This 
is merely a reflection of the fact that the charge distributions in the two layers will be 
relatively shifted, because the electrons intend to tunnel along the direction of the total 
magnetic field. If the electrons of the DQWS are in the symmetric state for B\\ = 0, where 
electron distributions in the two layers are directly on top of each other, the inclusion of an 

8 



in-plane field damages the original interlayer correlations and results in an increase in the 
interlayer Coulomb energy. In a pseudo-spin description |18[ of the DQWS, the tunneling 
behaves like a tumbling magnetic field, which twists the interlayer phase coherence of the 
symmetric states. In an attempt to minimize the total energy, a DQWS is forced to choose 
between the hopping energy and the Coulomb energy. The electrons can only take advantage 
of the hopping energy at the cost of increasing the interlayer Coulomb interaction. 

Putting Eq. ( PH) and Eq. (|i~2l) into Eq. (||), one obtains the desired equations of motion 
for Gn and G21 in Matsubara frequencies 

- -r £[^ qAq ' /2 £>n(q - q')<5ii(q', ico n ) + e^'^q - q'-)<5 21 (q'-, ico n )} 

+(iku n + /i)Gu(q, iu n ) + te- k2 ^ i e- u ^ kB/2 G 21 (q - k B , iu n ) = M q , , 
£[e^ Aq ' /2 t7 2 i(q' - q-)G n (q!, iu n ) + e^<™'~/ 2 f7n(q - q')G 21 (q'- ,iu n )] 

o q/ 

+(ihu n + /i)G 21 (q-, iu n ) + te- k B l o/ 4 e a ^ kB / 2 G 11 (q- + k B , iu n ) = 0, (13) 
where we have adopted the following notations 

^(q) = e^ a / 2 ^(q) 
G y (q) = e-^ 2 Gy(q) 
<p ij (q) >= e- iqa/2 < Plj (q) > 

q ± = q ± k B or q. (14) 

Eq. (^) and Eq. ( |14"D are intended to be applicable to different possible WC states. We 
need to set correct conditions for each situation. For the case of a perpendicular magnetic 
field, one has kg = 0, q = q, and a as the relative shift of the charge distributions between 
the layers. As mentioned earlier, there are two possible phases in the presence of an in- 
plane field, electrons can choose either to take advantage of the hopping energy at the cost 
of increased Coulomb energy, or to maintain good correlations for reducing the Coulomb 
energy at the cost of giving up the hopping energy. For the former case, one has = q±k B , 
and a = ksll^- For the latter case, one has q — q and can effectively set t = (see Eq. 



Eq. ([14]) can be rearranged into following compact matrix form, which is convenient for 
numerical evaluation, 

[(ihuj n + n) — D]G = HB (15) 

where 

B= [1,0,0,0,0, ] 

G= [G 11 (qi),G 21 (qr),G 11 (q 2 ),G 21 (q 2 -),...] (16) 

with {qi,q2,q3, •••} arranged in the order of increasing magnitude. The non-zero elements 
of the Hermitian matrix D are 

D 2i -i,2i = -t e -WV«2fe*B/2 
p 2 

L>2i-i,2j-i — —j-e u u\ c Li ~ Qj'J 

2 

n — Jftn^ij / 2 n* (n „-\ 
J^2i-i,2j — —j-e 3 <J 2 i{(\i — ) 

D 2i ,2i-i = -te-^o/i e ^ B /2 

D2i, 2j = ^-e^^Uui* ~ qj) (17) 
Eq. ( p~5| ) can be solved by diagonalizing the matrix D. If V k and uo k are respectively the 



fc-th eigenvector and eigenvalue of D, we obtain JTJ 



kmax 

<piMi) >= E v k (2i - i)v k (i) 

k 

< M<U) >= E V k (2i)V k (l) (18) 

k 

where k max is determined by fixing the chemical potential 

< Pll (0) >= u/2. (19) 

With the order parameters known, we can obtain the ground state energy per electron, 
which will be shown in the next section, from Eq. (§), 
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Ue-^[< p 21 (k B )> + < 

^ E{%")l<Mq")> 



Pi2(-ks) >] 

2 



z/eZ 



■[K(g) - V b (q) + K(?)cos(q ■ a)] | < p n (q) > | 2 }. 



(20) 



III. NUMERICAL RESULTS AND DISCUSSIONS 

We now discuss our numerical results for the different WC states. For this discussion, 
we compare the energy of the different WC states and then find the phase diagrams as the 
sample parameters are changed. The order parameters < Pij(q) > are obtained from a 
numerical analysis on Eq. ([T5|), in which well convergent results are obtained by keeping 16 
(24) shells in reciprocal-lattice vectors for triangular (square) WC states. In the following, 
we will first discuss the situation with a perpendicular field and then the situation with a 
tilted field. 

There are several possible configurations for a WC in a strong perpendicular field for 
different layer separations and hopping strength. At small layer separations, the electronic 
states of the different wells mix to form symmetric and anti-symmetric states. In the ideal 
case of d = 0, all electrons reside in the symmetric state for any value of hopping. A DQWS 
under this condition behaves as a single layer system. A WC in the symmetric state has a 
lattice constant a a such that {\f?>/2)a 2 = 2ti1 2 ./v. The order parameters at corresponding 
RLVs are 

< Pss (G) > = -£e-^+^ 2 cW +GH < C\ a C sP >^ 0, 

9 a/3 

<Pij(G) > = - <Pss(G) > for i,j = 1,2, (21) 

where C s(a) = (1/ V / 2)(C 1 ± C 2 ). The above expression also shows that the charge dis- 
tributions for the two layers are directly on top of each other, i.e., a = 0. At finite, but 
small, layer separations, the WC state of a DQWS is essentially the above one-component 
triangular (OCT) lattice. For large enough layer separations, the symmetric state is no 
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longer energetically favored, as the system begins to seek a state where electrons within the 
same layer are more strongly correlated than electrons in different layers. In the large d 
limit, a DQWS becomes two independent single layers. Electrons in each well form their 
own triangular WC. These two-component triangular (TCT) lattices have a lattice constant 
(v^3/2)a„ = 2itIq/(v/2), larger by a factor of \[2 than that of the OCT lattice discussed 
above. To minimize static interlayer Coulomb energy, the two WCs are relatively shifted so 
that the lattice sites of one WC are at the centers of the triangles of the other WC lattice, 
i.e., a = (l/3)(a + b Q ), where a Q and b Q are the primary lattice vectors. A shifted two- 
component square (TCS) lattice WC state with a = (l/2)(a Q + b D ), where a G and b Q are 
the primary lattice constant of the square WC, has lower interlayer Coulomb energy than 
the above TCT WC state, since the lateral distance between a electron in one layer and its 
nearest electron in the other layer in a shifted TCS lattice structure is larger than that in 
a shifted TCT lattice structure. For an intermediate range of layer separations, this TCS 
lattice structure may become the ground state of a DQWS. 

In Fig. ([]]), we show the energies of the three different WC structures discussed above 
as functions of layer separation d for different value of hopping t. The lowest energy states 
are the OCT WC at small d, the TCS WC at intermediate d, and the TCT WC at large d. 
The range of layer separations at which the TCS WC exists as the ground state shrinks, and 
finally disappears, when t is increased, as a result of the expansion of the OCT WC phase. 
The important conclusion from Fig. (|l|) is that for weak interlayer hopping, two structural 
phase transition should be expected when the layer separation is increased: first from a 
OCT WC to a TCS WC, then from a TCS WC to a TCT WC. 

Next, we consider the influence from an in-plane magnetic field. Our discussion will 
concentrate on the situation where the system was in the symmetric state prior to the 
application of the in-plane field, since the effects of an in-plane field are the largest in the 
symmetric state due to the existence of strong inter-layer coherence. As a result of the 
competition between hopping energy and Coulomb energy, a DQWS in a tilted field can be 
in either of the two different ground states depending on the value of the in-plane field. One 
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is the symmetric state WC (SSWC) described by Eq. (|21|), where a = and the electrons 
feel no effect of the hopping and B» has no effect. The other, which is generally relevant at 
a small in-plane field, is the state in which electrons form linear combinations of the states 
in the two wells that are displaced by the total magnetic field. We call this state a twisted 
symmetric state. Defining 

C-s(a)a — — [Cla ± CW-fcg]; (22) 

the order parameters for a twisted symmetric state WC (TSSWC) are 

< e „(G) >= - E e- lG * (a+f3)/ %, a+Gy < C\ a C sP 0. (23) 

9 a/3 

The order parameters in layer representation < Pij(G) > for this TSSWC can be obtained 
from the above expression by making use of the operator relation of Eq. (p2|), 

< pn(G) >= - < g ss (G) > 
<p 22 (G) >= ^ G * kBl *° < Qss(G) > 

< Pl2 (G - k s ) >= l -J G * kBl ° < Qss (G) > 

< p 21 (G + k B ) >= i e * G ^° < Qss (G) > . (24) 

From the above expression, one can see the obvious effects of the in-plane field: it shifts 
the charge distributions of the two layers relatively by an amount a = ksl 2 ^ = d(B\\/B±)5t; 
The positions of the non-zero interlayer order parameters < Pi2(2i)(q) > m the reciprocal 
vector space are shifted from G to G = G ± kg. The result is an increase in interlayer 
Coulomb energy. At small values of in-plane field, the increase in the Coulomb energy 
is small and can be compensated by the hopping energy, so the TSSWC is favored over 
the SSWC. As the in-plane field increases, the cost in Coulomb energy increases while the 
hopping energy decreases. When the gain in the hopping energy can no longer compensate 
the cost in Coulomb energy for a strong enough in-plane field, the DQWS undergoes a 
transition from the TSSWC to the SSWC, where the total energy of the system is lowered 
by giving up the hopping energy to restore the original interlayer coherence. In Fig. (§), we 
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show the energies of the TSSWC and SSWC as functions of the in-plane field for different 
values of hopping energy. It is clear, from the figure, that the TSSWC is the energetically 
favored ground state at small in-plane field while the SSWC is the energetically favored 
ground state at large in-plane fields. The critical value of the in-plane field for this phase 
transition as a function of the hopping energy is shown in Fig. (|3|), larger critical values of 
B\\/B± for larger values of hopping t. 

Since a WC breaks the rotational symmetry of the system, the properties of a DQWS, in 
principle, depend on the angle between the in-plane field and the crystal axes. As the WC 
is pinned by the presence of weak disorder, the angle can be changed by simply sweeping 
the direction of the in-plane field. This provides a potential opportunity to probe the 
orientational order of a WC and to find a unambiguous signature of the existence of a WC. 
In Fig. (|), we show, for t = 0.02, t = 0.08, and t = 0.4, the energies of the TSSWC for 
both the cases where the in-plane field is perpendicular to (the solid-lines), or parallel with 
(the dot-lines), one of the crystal axes. We can see that the differences are small. The 
change in the value of critical in-plane field from the different orientations of the field is 
practically indistinguishable. This is mainly because that the phase transition occurs at 
an in-plane field where |a| is small compared to the lattice constant for reasonable sample 
parameters. However, we should not rule out the possibility that some other quantities 
may have a measurable dependence on the anisotropy of the WC. For example, the energy 
difference of the TSSWC for the two different orientations of the in-plane field right before 
the phase transition is As ~ 7 x 10~ 4 (e 2 /eZ o ) for t = 0.08(e 2 /eZ o ), which is on the order of 
70mK for e 2 /el ~ 100K. IF the same order of anisotropy exists in the low-lying phonon 
modes (see the next section), we may have a measurable difference in the anisotropy of the 
specific heat at low enough temperatures, say, T < 70mK. 
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IV. CRYSTAL PROPERTIES IN A PARALLEL MAGNETIC FIELD 



In this section, we discuss some consequences of immersing the double well system in a 
parallel magnetic field, when the electron system is in the TSSWC - i.e., a lattice of particles 
in symmetric states of the two wells, with the position of the single-particle orbits in one 
well displaced in the direction of the total applied magnetic field. As was discussed in the 
previous section, and as shown in Fig. 2, the energy per particle in the electron lattice 
depends (weakly) on the orientation of the crystal axes relative to the parallel magnetic 
field. One may think of this effect as an energy cost for having the bond angles of the lattice 
deviate from some preferred direction. While the energy cost per particle may be small, long 
wavelength fluctuations of the lattice - where many bond angles deviate from the preferred 
direction - will be strongly affected by this energy cost. This will have important effects on 
the long wavelength collective modes of the system, as well as the state of the crystal in the 
presence of a slowly varying random potential (which typically arises in real heterostructure 
environments.) 

To model the bond-angle energy, we use a continuum elasticity theory approach. The 
energy of a two-dimensional lattice deformed from a perfect crystal configuration by a dis- 



placement field u(r) may be written as |20 



Z J CLc\ 



L kk)i 

where /i and A are Lame coefficients, a® is the lattice constant, Uij = d{Uj + djUi is the strain 
tensor, and repeated indices are summed over. The bond angle field may be written in terms 
of the displacements in the form 

6{r) = -(d x u y -d y u x ), 

so that it is natural to write the energy of a configuration in the presence of parallel magnetic 
field in the form E = E + Eg, with 



E e = i eo J d 2 r9(f) 
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where eo is a phenomenological parameter describing the energy cost to misalign the crystal 



with the parallel magnetic field pl|. 



This model has been studied in the context of a two-dimensional crystal adsorbed on a 
periodic substrate , where it was noted that the extra term tends to increase the stability 
of the crystal with respect to finite temperature. It is also interesting to note that such a 
bond-angle term increases the stability of the crystal with respect to quenched disorder. To 
see this, consider a weak random potential acting on the WC, which has a long orientational 
correlation length £ (or even quasilong-range orientational order ||19|| ). The energy cost to 
misalign a correlated region with the magnetic field scales as £ 2 due to Eg, whereas the 
energy gained by aligning with the random potential scales only as V^c ~ where N c is 
the average number of electrons in a correlated region. Thus, for weak enough disorder, 
where £ is large, one expects the correlated regions to align with the preferred orientational 
axis. 

An interesting possible method to demonstrate this would be to anneal a WC in the 
double well system in the presence of a parallel magnetic field. It has been shown recently 
|19| that a WC in the presence of a slowly varying random potential freezes into a state 
with at best power-law (i.e., quasilong-range) orientational order. By cooling the system in 
the presence of the parallel field, a preferred orientational axis is picked out for the crystal, 
leading to the possibility that one could obtain a system with true long-range orientational 
order. Observing that phenomenon will be possible if the orientational correlation length is 
large enough at the freezing transition that the Eg term overcomes thermal fluctuations in 
the orientation of a correlated region of the lattice; i.e., e N c > T M , where T M is the melting 
temperature of the crystal. Experimentally, one could probe this effect by measuring the 
depinning electric field of the lattice, which has been shown to be sensitive to orientational 
correlations of the crystal one expects to see a decrease in the pinning field if the 



orientational correlations are increased. The possibility of creating a WC with long-range 
order using parallel magnetic fields in a double well system is currently under investigation 
23|. 
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It is also interesting to consider the effect of the bond-bending term on the collective 
mode spectrum of the WC. It is well known that in a perpendicular magnetic field, the 
WC supports a phonon mode dispersing as uj(q) oc q 3 ^ 2 , and there have been attempts to 
measure this directly using rf absorption |J. Since the bond-angle term Eg represents a 
restoring force on long-wavelength fluctuations, we expect the phonon energy to change in 
the present condition. For this purpose, we write down the appropriate form for the energy 
of a crystal deformation in terms of a dynamical matrix that yields the energy E = Eq + Eq 
in the long wavelength limit. The low frequency collective mode frequencies are given by 
X a /u c , where X a are the eigenvalues of the matrix a y D, o y is the Pauli spin matrix, and 
D = D° + D e is the total dynamical matrix [2f|. D°, the dynamical matrix in the absence 
of a parallel magnetic field, has been shown to be 



q 

where b and A^ a ^ are constants. The dynamical matrix associated with Eq is easily obtained, 



q a/3 



D e Jo 2 D e = Jo 2 D s =D e = _ Jo 

4m v ' yv Am y v 4m 

By noticing the similarity between D e and the second term in D°, it is easy to show that 
under the present condition, the low- lying phone mode still disperses like uo(q) = Cq 3 ^ 2 , but 
the coefficient C is increased, with a new contribution from the dynamic matrix D e . 



V. CONCLUSION 

Working in the Hartree-Fock approximation, we calculated the energy of different Wigner 
states of the two-dimensional electron gas for a double quantum well system in a strong 
magnetic field. We found the phase diagram for the evolution of the WC states in a strong 
perpendicular magnetic field when layer separation and hopping are changed. In the ab- 
sence of interlayer hopping, the ground state at small layer separations is a one-component 
triangular lattice Wigner state which possesses interlayer coherence. As the layer separation 
is increased, the ground state first undergoes a transition to a two-component square lattice 
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Wigner state, and then undergoes another transition at an even larger layer separation to a 
two- component triangular lattice Wigner state. The range of the layer separations at which 
the two- component square lattice occurs as the ground state shrinks, and eventually disap- 
pears, as the interlayer hopping is increased. We also studied the in-plane magnetic field 
induced phase transition in the Wigner state, which has so far only been studied experi- 
mentally for v = 1 quantum Hall state. We calculated the critical value of the in-plane field 
for the transition. We find that the anisotropy of the Wigner state, i.e.,, the orientation 
of the crystal with respect to the direction of the in-plane magnetic field, has a negligible 
effect on the value of the critical in-plane magnetic field for small filling factors. The effect 
of this anisotropy on the low-lying phonon energy is discussed. A possible experimental 
arrangement for observing the in-plane field enhancement of the orientational order in the 
crystal in the presence of a weak disorder potential is also discussed. 
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FIGURES 

FIG. 1. Energy per particle for a one-component triangular lattice (OCT), a two-component 
square lattice (TCS), and a two-component triangular lattice (TCT) Wigner crystal as functions 
of layer separation d for v = 1/4 at different values of hopping t. (a) t=0; (b) t = 0.01(e 2 /d o ); (c) 
t = 0.02(e 2 /eZ o ). 

FIG. 2. Energies per particle of a twisted symmetric state Wigner crystal (TSSWC) and a 
symmetric state Wigner crystal (SSWC) as functions of in-plane field at v = 1/4 and d = l Q for 
different value of hopping. The dash-dot line is for the SSWC. The solid lines are for the TSSWC 
with the in-plane filed perpendicular to one crystal axis. The dot-lines at t = 0.02, 0.08, 0.4(e 2 /e/ o ) 
are for the TSSWC with the in-plane field parallel with one crystal axis. 

FIG. 3. The critical value of B\\/B± as a function of the hopping for v = 1/4 and d = l Q . 
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